Probing the superconducting gap symmetry of PrOs4Sbi2: A penetration depth study 
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We report measurements of the magnetic penetration depth A in single crystals of PrOs4Sbi2 
down to ~0.1 K using a tunnel-diode based, self-inductive technique at 21 MHz, with the ac field 
applied along the a, b and c directions. In all three field orientations the penetration depth and 
superfluid density tend to follow a power law, A ~ T^jp^ ~ at low temperatures, ps for 
various gap functions have been calculated, and data is best fit by the ^He A-phase-like gap with 
multidomains, each having two point nodes along a cube axis, and parameter Ao(0)/A;b Tc = 2.6. 
This suggests that PrOs4Sbi2 is a strong-coupling superconductor with two point nodes on the 
Fermi surface. We also confirm the observation of the double transitions at 1.75 K and 1.85 K seen 
in other measurements. 



The recent discovery of the Heavy Fermion (HF) skut- 
terudite superconductor (SC) PrOs4Sbi2 P, has at- 
tracted much interest due to its differences with the 
other unconventional SC, and in particular, the HFSC. 
PrOs4Sbi2 has a nonmagnetic ground state of locahzed 
/ electrons in the crystalline electric field, hence its HF 
behavior, and consequently the origin of its supercon- 
ductivity, might be attributed to the interaction between 
the electric quadrupolar moments of Pr'^"'" and the con- 
duction electrons. It is thus a candidate for the first 
SC mediated by quadrupolar fluctuations, i.e. by nei- 
ther electron-phonon nor, as with other HFSC, magnetic 
interactions. 

Recent experiments on PrOs4Sbi2 give conflicting ev- 
idence to the nature of the SC gap. Muon-spin rotation 
(/iSR) measurements revealed a low-temperature expo- 
nential behavior, suggesting isotropic pairing (either s or 
p-wave) 0. Scanning tunneling spectroscopy measure- 
ments also measured a density of states (DOS) with no 
low-energy excitations with a well-developed SC gap over 
a large part of the Fermi Surface (FS) 4]. The absence of 
a Hebel-Slichter peak and the non-T'^ behavior of 1/Ti in 
nuclear quadrupolar resonance (NQR) experiments sug- 
gest that PrOs4Sbi2 has a full gap or point nodes, but 
not line nodes, at zero fleld [^|^. If PrOs4Sbi2 has an 
isotropic gap, then it is unique among HFSC, suggesting 
the possibility of (a) an important difference in super- 
conducting properties between HFSC with magnetic and 
non-magnetic /-ion ground states, and (b) a correlation 
between pairing symmetry (isotropic or nodal gap) and 
mechanism (quadrupolar or magnetic fluctuations) of su- 
perconductivity 3] . 

Unlike the /iSR and NQR results, angle-dependent 
thermal conductivity measurements Q revealed two dis- 
tinct SC phases with different symmetries, a phase tran- 
sition between them, and presence of point nodes. In the 
high-field phase four point nodes ([100] and [010] direc- 
tions) have been observed, whereas there are only two 



point nodes ([010] directions only) in the low-field phase. 
Specific-heat data Q also show a low-temperature power 
law behavior, suggesting the presence of nodes. An- 
other recent /xSR experiment revealed the spontaneous 
appearance of static internal magnetic fields below Tc, 
providing evidence that the SC state is a time-reversal- 
symmetry-breaking (TRSB) state 0. 

In this paper, we present high-precision measurements 
of penetration depths A(r) of PrOs4Sbi2 at tempera- 
tures down to 0.1 K. The ac field was applied along all 
three crystallographic axes. In all three field orientations 
both A(r) and superfluid density Ps{T) tends to follow 
a quadratic power law, suggesting that the SC gap has 
nodes on the FS. ps for various gap functions has been 
calculated, and data are best fit by the He-3 A-phase-like 
gap, with two point nodes in the [010] directions. Our 
data thus puts PrOs4Sbi2 in line with other HFSC, in 
that they all have nodes on the FS, despite the proposed 
non-magnetic nature of the mechanism of its supercon- 
ductivity. 

Details of sample growth and characterization are de- 
scribed in Ref.'^. Measurements were performed utilizing 
a 21-MHz tunnel diode oscillator (1^] with a noise level 
of 2 parts in 10^ and low drift. The magnitude of the ac 
field was estimated to be less than 5 mOe. The cryostat 
was surrounded by a bilayer Mumetal shield that reduced 
the dc field to less than 1 mOe. The sample was aligned 
inside the probing coil in all three crystallographic direc- 
tions. The sample was mounted, using a small amount of 
GE varnish, on a single crystal sapphire rod. The other 
end of the rod was thermally connected to the mixing 
chamber of an Oxford Kelvinox 25 dilution refrigerator. 
The sample temperature is monitored using a calibrated 
Ru02 resistor at low temperatures {Ttase - 1-8 K), and 
a calibrated Cernox thermometer at higher temperatures 
(1.3 K - 2.5 K). 

The deviation AA,(T) = Xi{T) - A,(0.1 K) {i = a, 
b, c) is proportional to the change in resonant frequency 
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A/(T), with the proportionahty factor G dependent on 
sample and coil geometries. The subscript i denotes the 
direction of the applied magnetic field. For our H//a 
data, we determine G from a single-crystal sample of pure 
Al by fitting the Al data to extreme nonlocal expressions 
and then adjusting for relative sample dimensions |lll| . 
Testing this approach on a single crystal of Pb, we found 
good agreement with conventional BCS expressions. The 
value of Ga obtained this way has an error of ±10 %, 
since our sample has a rectangular basal area instead of 
square ^^]. To obtain Gf, and Gc we make use of the 
cubic symmetry of the crystal and assume that the total 
change in penetration depth from the three orientations 
are equal, i.e. AAa(Tc) = AAf,(Tc) = AAc(Tc). From 
this equality, and the value of Ga^ we can calculate Gb 
and Gc- 

Fig. n shows AAi(r) as functions of temperature. All 
three curves vary strongly at low temperatures, incon- 
sistent with exponential behavior expected for isotropic 
s-wave superconductors. On the other hand, the varia- 
tion is not linear, but has an obvious upward curvature, 
unlike the low-temperature behavior expected for pure 
(f-wave superconductors. A fit of the low temperature 
data (up to 0.55 K w 0.3 Tc) to a variable power law, 
AA(T) = ^ -I- ST" yields n = 1.9 ± 0.1 for H / / a, b, 
and n — 2.0 ± 0.1 for H / / c. The uncertainty in n is not 
a consequence of the uncertainty in G, but rather of the 
somewhat uneven faces of the crystal and the range of fit. 
Within the uncertainty in G the three curves are consis- 
tent with one another. There is also a small upturn near 
0.62 K, which might distort the low-temperature fit and 
cause the power law to deviate from . The NQR spin- 



FIG. 2: (a) Low-temperature superfluid density ps[T) = 
[X^{0)/X^{T)] calculated from AA(r) data in Fig. El for all 
three field orientations. (Q) H// a, (A) H//b, (x) H//c. 
Using Ao(0)/fcfl Tc = 2.6, the solid lines are the calculated ef- 
fective superfluid density pl^-^ corresponding to gaps II, III, A 
and B. The dotted lines correspond to ps{H / /b) and Ps{H / /a, 
c) for gap B. (b) The same calculated curves over the entire 
temperature range. Only the Ps{T) data for H j ja are shown. 



lattice relaxation rate also changes around this temper- 
ature, however, the origin is not clear at present [sl Il3|. 
A non-unitary state has the unique feature that spin- 
up and spin-down Cooper pairs have different excitation 
gaps [3. If the SC state in PrOs4Sbi2 is a TRSB state, 
then this upturn may be due to the contribution from 
the smaller gap It is interesting to note that a fit 
of the AA(r) from 0.6 K to 1.1 K, to the same vari- 
able power law, gives an exponent of about 3, consistent 
with n « 4 obtained for specific-heat data over the same 
temperature range 0. 

Using the value of A(0) — 3440 A from ^SR measure- 
ments we calculated the superfluid density ps from 
our data. We follow the procedure in Ref . O to compute 
the experimental superfluid density, using the fit to 
estimate the small difference between A(0) and A(0.1 K). 
Fig. shows ps ( for all three field orientations at low 
temperatures. In each case, a fit of PsiT) to a variable 
power law, Ps{T) = 1 - aT" also yields n « 2, from 
0.1 K (~ 0.05 Tc) to 0.55 K. Once again this suggests 
the presence of low- lying excitations, incompatible with 
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FIG. 3: Polar plots of gap functions (II), (III), (A) and (B). 

an isotropic SC gap. 

Several theoretical proposals have been put forward to 
understand the two SC phases 0, 0, Q|. To explain 
the behavior of the angle-dependent jinagneto-thermal- 
conductivity results 0, Maki et al. [l^l proposed three 
possible SC gap functions for PrOs4Sbi2. In particu- 
lar, for the low-field (L) phase, two gap functions were 
proposed: (II) /(k)=l - fc^ - fc^ having four nodes and 
(III) /(k)=l - ky, with two point nodes. The gap 
function is A(k) = Ao/(k), with the form factor /(k) 
normalized to unity and Aq the temperature- dependent 
maximum gap value. As we will see, both functions 
lead to a linear temperature dependence for the super- 
fluid density. Consequently, we consider two further 
gap functions: (A) d(k)=y ka (a line-node gap), and 
(B) d{k)=y {ka± ifcc), for which A(k)= Aq |<i(k)| . 
Gap B has two point nodes along the [010] directions 
and a gap dispersion identical to the superfluid ^He A- 
phase, A(k) = Aq \ka ± ikb\ = Aosin6'; hence, they give 
identical temperature dependences of ps- Polar plots 
of these gap functions are shown in Fig. We have 
assumed the gap maximum Ao(T') to have the form 

Ao(T) = S,ckTctanh{-^^a{^){^ - 1)} [ItI, where 
Ssc = Ao(0)/A;b Tc is the only variable parameter, Tc = 
1.85 K, o = 2/3, and the specific heat jump AC/C = 3 
is an experimentally obtained value [18l |. 

A problem arises immediately with point nodes. If 
there are only two point nodes in the [010] directions, 
breaking cubic symmetry, then one would expect AXb to 
tend toward an exponential temperature dependence at 
low temperatures. We show this in Fig. [3 where we have 
calculated the superfluid density for gap B for fields along 
[010] and either [100] or [001]. A measurement along 
[010] would indeed give exponential behavior while mea- 
surements in orthogonal directions give a strong temper- 
ature dependence. However, our experimental data show 
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FIG. 4: Quasiparticle DOS for the gap functions II (full 
line). III (short-dashed line), A (short-dash-dotted line) and 
B (short-dotted line). 

otherwise - there is an almost identical superfluid re- 
sponse in all three field orientations. While it is possible 
that, in the absence of an external agent, the sample will 
randomly choose one, and only one, axis along which to 
locate the nodes each time it becomes superconducting, 
it is much more likely to develop a domain structure. One 
possibility is that the SC order parameter and strain are 
coupled [iq . a situation similar to magnetostriction in 
ferromagnets. Such a situation also arises in Cromium 
(Cr) [ij], where the coupling between the spin-density- 
wave (SDW) and strain wave causes the wave vector Q of 
the modulation to point along any {001} direction in the 
bcc Cr lattice. In bulk Cr all three possible orientations 
occur with equal probability. A Cr single crystal thus 
has multiple domains, with each domain corresponding 
to one of three possible Q^^, Qy and regions - 
the "poly-Q" state. Tensile stress in applied along one 
cube axis while cooling through Tpf, however, produces a 
"single-Q" state with all domains having their Q's point- 
ing in the same direction, along the stress axis [21| . Evi- 
dence that a domain with only two point nodes can form 
in PrOs4Sbi2 was reported in Ref. where nodes were 
seen only along a single [010] axis. We suggest that the 
experimental setup of Ref. 7 may have strained the sam- 
ple to produce a single domain, analogous to the single-Q 
state of Cr. 

Assuming the existence of domains, we plot an effec- 
tive superfluid density p'^^^ by taking the average of [100] , 
[010] and [001]-superfluid densities, with equal weight 
from each component. The superfluid densities in dif- 
ferent directions are evaluated using the expression 

P.w/^) = i-^E(^.^ + ^?)^ (1) 
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where x, y, z = any permutation of a, h, c. Nifi) is the 
quasiparticle DOS at the FS, / = [exp(Ej^/fcsr)+l]-i is 
the Fermi function, and = [£2(k) + A(k, r)2]i/2 jg the 
quasiparticle energy. The component superfluid densities 
for Gap B are shown as two dotted Hues, and p*^/^ as a 
sohd hne, in Fig. [3 Clearly, the agreement between data 
and pI^^ is very good. We chose the strong-coupling 
value 5sc = 2.6 here, taken from Ref. 0. Using = 2.1 
from Ref. Ogives a worse fit. For the other gap functions, 
we also calculated p%^^ (shown in Fig.|21) and ps{H//a, b, 
c) (not shown here) - all of them give linear temperature 
dependences and fall far from the experimental data. The 
effective quasiparticle DOS for all four gap functions are 
also shown in Fig. 01 Our data therefore suggest that 
PrOs4Sbi2 is a strong-coupling unconventional SC. The 
superfluid data is best fit with a ^He A-phase-like gap, 
with two point nodes on its FS. Both the field-direction- 
independence of the superfluid data, and the nice fit of 
the data to Ps^-^ , strongly suggest that PrOs4Sbi2 has 
multidomains. Note that thou gh Gap B is a unitary gap, 
our low-T data can also be fit [22| by the two-point-node 
L-phase non-unitary gap proposed by Ichioka et al. in 
Ref. 0, with similar DOS structure. Hence our data does 
not contradict the non- unitary result of Aoki et al. j^. 

It is already apparent in Fig. that the data devi- 
ate from the proposed gap function above 0.6 K. This is 
even clearer in Fig.l^b, which shows Ps{T) for H / / a from 
0.1 K to Tc- The other two field orientations give almost 
identical temperature dependence. None of the four gap 
functions fits the data over this larger temperature range. 
This could be due to the opening up of the smaller gap 
caused by the non-unitarity of the SC state, as mentioned 
earlier. Also, strong changes in the mass renormalization 
in different sheets have been found in de Haas van Alphen 
experiments 0. These changes may cause the distribu- 
tion of values of the SC gap measured in tunneling mea- 
surements 0, and strengthens both the idea that the 
mass renormalization and superconductivity are of the 
same origin, i.e. that the quadrupolar fluctuations favor 
SC correlations, as well as the possible multiband char- 
acter of superconductivity in this material jj,]. Hence a 
multiband analysis, similar to those performed on MgB2, 
might be required to fit the superfluid data over the en- 
tire temperature range. 

Finally, we turn to the region near Tc- Fig. [S] shows 
AA(T') near Tc for all three field orientations. Three 
features can be seen: the onset of superconductivity at 
1.88 K, a strong but broad shoulder near 1.8 K, and fi- 
nally a weak shoulder near 1.74 K (observable even in the 
H 1 1 c data). In another sample from the same batch, 
only the first and third features were observed. The 
1.88 K and 1.74 K features confirm the two superconduct- 
ing transitions seen in the specific- heat measurement , 
and suggested by angle-dependent thermal conductivity 
measurements J\. The origin of the 1.8 K shoulder is 
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FIG. 5: AA(r) for all three field orientations near Tc- The 
curves have been offset for clarity. The arrows indicate the 
second superconducting transition at 1.74 K. 



unknown. In the Ps{T) plot, ps already approaches zero 
near 1.7 K. So these features were not discernible there. 
Also, we did not see any anomaly around T* = 2.3 K 
that was observed in Ref. 0. 

In conclusion, we report measurements of the magnetic 
penetration depth A in single crystals of PrOs4Sbi2 down 
to 0.1 K using a tunnel-diode based, self-inductive tech- 
nique at 21 MHz, with the ac field applied along the a, 
b and c directions. In all three field orientations A and 
superfiuid density ps tend to follow a quadratic power 
law. We have calculated Ps for various gap functions, 
finding that the data are best fit by the '^He A-phase-like 
gap function with two point nodes on the FS. We also 
observe the double transitions near 1.75 K and 1.85 K 
seen in other measurements. 
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